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We study the zero temperature phase diagram of hardcore bosons on the dual of the bowtie lattice. 
Two types of striped diagonal long-range order (striped order I and striped order II) are discussed. A 
state with type-II striped order and superfluidity is found, even without nearest-neighbor repulsion. 
The emergence of such a state is due to the inhomogeneity and the anisotropy of the lattice structure. 
However, neither the translational symmetry nor the symmetry between sublattices of the original 
lattice is broken. In this paper, we restrict a 'solid state' of lattice bosons as a diagonal long-range 
ordered state breaking either the translational symmetry of the original lattice or the symmetry of 
different sublattices. We thus name such a phase a striped superfluid phase (SSF). In the presence 
of a nearest-neighbor repulsion, we find two striped charge density wave phases(SCDW I and II) 
with boson density p = 1/2 (with striped order I) and p = 2/3 (with striped order II) respectively, 
when the hopping amplitude is small compared with the repulsion. The SCDW I state is a solid, 
in which the translational symmetry of the original lattice is broken. We observe a rather special 
first-order phase transition showing an interesting multi-loop hysteresis phenomenon between the 
two SCDW phases when the hopping term is small enough. This can be accounted for by the special 
degeneracy of the ground states near the classical limit. The SSF re-appears outside the two SCDW 
phases. The transition between the SCDW I and SSF phases is first order, while the transition 
between SCDW II and SSF phases is continuous. We find that the superfluid stiffness is anisotropic 
in the SSF states with and without repulsion. In the SSF with repulsion, the superfluid stiffness is 
subject to different types of anisotropy in the region near half filling and above 2/3-filling. 
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The supersolid, a novel quantum state with diagonal 
simultaneous long-range order(DLRO) and ofF-diagonal 
long-range order(ODLRO), was introduced nearly half a 
century ago- . This exotic quantum state has attracted 
considerable research interests in recent years. Owing 
to the fast development of laser cooling techniques, var- 
ious optical lattices can be realized^, which makes the 
investigation of supersolidity of bosons on discrete lat- 
tices more realistic. For softcore bosons, the supersolid 
phase emerges by the " defect condensation" mechanism, 
where doped bosons(holes) act as interstitials (vacancies) 
in the crystal3^-S. However, for hardcore bosons with 
nearest-neighbor repulsion, no supersolidity was found 
on the square^i^, honeycomb^, kagome^, star latticeiS 
and Shastry-Sutherlancl^ lattices, due to the instability 
of such a phase^'^ on these lattices, which leads to phase 
separation into a pure solid and a superfluid for all val- 
ues of the interaction strength. The situation changes in 
the presence of geometric frustrationi^r— . Supersolidity 
was found for hardcore bosons on the triangular lattice, 
where the extensive degeneracy of the classical ground 
states at half filling is lifted by quantum fluctuations and 
the ground state attains DLRO and ODLRO simultane- 
ously, thus forming a supersolid. Moreover, it was sug- 
gested that frustration induced by next-nearest neighbor 
interactions^^ can also lead to supersolidity. These find- 
ings show a different route to supersolidity, which is based 
on an order-by-disorder mechanism, by which a quantum 
system avoids classical frustration. However, all these in- 
vestigations considered simple homogeneous lattices. It 
is natural to investigate the behavior of bosons also on 



inhomogeneous lattices. 

Let us consider hardcore bosons with nearest-neighbor 
repulsion on the dual of the bowtie lattice, which is 
anisotropic in two directions and consists of two types of 
sites with different coordination number, see Fig. [TJ The 
corresponding Bravais lattice is a rectangular lattice. It 
has a basis of 3 sites A, B, and C with the primitive 
vectors ai , sl2 as shown in Fig. [TJ The Hamiltonian is 

where a] {ai) creates (annihilates) a boson at site i, t is 
the nearest-neighbor hopping amplitude, V the nearest- 
neighbor repulsion, /x the chemical potential, and Ui = 
or 1. 

The model can be mapped onto the spin-1/2 XXZ 
model on the same lattice with Sz{i) — rii — 1/2 in the 
usual way, 

(2) 

where J = 2t is the in-plane exchange, = V is the ex- 
change in the z direction, and h{i) is a staggered external 
magnetic field given by h{i) — ^i — 2V for sites with coor- 
dination number 4 and /i — SV/2, otherwise. This marks 
the special character of the model. The solid state of the 
bosons is equivalent to magnetic order in the z direction. 

In the classical limit (i = 0), at zero temperature, there 
are only two phases if the nearest-neighbor repulsion (or 
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the exchange in z direction) V is absent. The lattice is 
empty when /i < 0, full if /i > 0. 

With the repulsion V present, there exists four phases. 
For /i/y < 0, the lattice is empty. For y^jV > 4, the lat- 
tice is full. In the region < n/V < 2, the lattice is half 
filled. The model shows a solid ordering (antiferromag- 
netic order in z direction) which breaks the translational 
symmetry of the original lattice, as shown in Fig. [Ija). 
Two bases half filled with three bosons form a new basis 
of a sohd at a wavevector (tt/S, 0), showing the character 
of a striped sohd or a SCDW^^. We refer to this order 
as striped order I. In the region 2 < fi/V < 4 the lattice 
is 2/3 filled, showing a DLRO again. Two of three sites 
in a basis are filled and form the basis of a phase with 
a wavevector (27r/3,0), as shown in Fig. [Ub). However, 
this ordering does not break the translational symmetries 
of the original Br avals lattice, nor the symmetry of A and 
B sublattices. (The spontaneously breaking of the sym- 
metry between two sublattices has been reported in @) . 
Thus, we would not say the bosons form a solid. In this 
paper, we restrict a 'solid state' of lattice bosons as a di- 
agonal long-range ordered state breaking either the trans- 
lational symmetry of the original lattice or the symmetry 
of different sublattices. The order is again of a striped or 
a SCDW type, to which we refer as striped order II. With 
quantum hopping present, SCDW phases were found in 
models with next-nearest neighbor repulsion^^ or plaque- 
tte interactions^^. For hardcore bosons on the dual of the 
bowtie lattice, without the next-nearest neighbor repul- 
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FIG. 1. (Color online). The dual of the bowtie lattice con- 
sisting of ni X 712 = 4 X 6 bases with periodic boundary con- 
ditions. The basis is shown in the top, where A, B and C 
label the three sublattices. The A and B sublattices are re- 
lated by symmetry, and ai and a2 are primitive vectors with 
length 3 and 2 respectively. The numbers arranged in horizon- 
tal (vertical) direction are the x{y) coordinates of the lattice 
sites. The linear size of the system is Lx = Ly = 12, with 
N = Lx X Ly /2 — 72 sites. The coordination number of sites 
connected by horizontal bonds is 3, and that of sites with 
an X coordinate satisfying mod (a;, 3) = is 4. (a) and (b) 
show the p = 1/2, 2/3 SCDW phases respectively, where black 
circles represent bosons and white ones stand for vacant sites. 




FIG. 2. (Color online). The possible boson (black dot) dis- 
tributions on the 24 X 4 (ri;, — 8) dual of the bowtie lattice at 
t = 0,fj. = 2V with density 12/24, 13/24, 14/24, 15/24, 16/24 
(from bottom to top, n from 4 to 0) respectively. 

sion or plaquette interactions, we shall show that striped 
phases also emerge. 

At fi/V = 2, a, special degeneracy of the ground states 
appears, as illustrated in Fig. [5] Since breaking one of the 
two stripe ordered patterns along the y-direction costs 
energy, no interface can be formed along the horizontal 
direction at zero temperature. The lattice can be divided 
into blocks along the x-direction without creating any 
such interfaces. The length of a half filled block in the 
x-direction is 6, that of the 2/3 filled one is 3. For a 
system with linear size L = Srif, in the x-direction (rib 
is the number of bases along the x-direction, which is 
restricted to be even in this work), let n and Ub — 2n be 
the number of blocks with particle density p = 1/2 or 2/3 
respectively, the energy density (per site) of the system 
is £■ = {-3nn+{nb-2n){-2fj, + V))/(3nb) ^ -V, which 
is independent on n. Thus, there are nb/2 + 1 possible 
boson densities: p — {2nb — n)/{3nb), n = 0, 1, • • • , nb/2. 
For each filling, there is a total of Nd{p) — nb{nb — n — 
l)!/(nfc — 2n)!/n! degenerate states. In other words, there 
are nb/2 + 1 phases coexisting. In the thermodynamic 
limit, the system can have any density between 1/2 and 
2/3. Although the zero temperature entropy per site is 
still in the thermodynamic limit, the long-range order 
is broken in these phases, except for the half filling and 
the 2/3 filling phases. In the language of the spin-1/2 
XXZ model, the long-range magnetic order is destroyed 
by the competition between the staggered field and the 
exchange in the z direction, instead of the geometrical 
frustration. This is different from the frustrated ground 
state on the triangular lattice at half filling. 

Considering these interesting properties of the model 
introduced by the inhomogeneity and the anisotropy of 
the dual of the bowtie lattice, it is desirable to explore 
quantum properties of the model. To this end, we per- 
formed extensive simulations by using the stochastic se- 
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FIG. 3. (Color online). Phase diagram at the limit V = 
0. The upper inset shows density wave: particle density per 
column as a function of the x coordinate of the column, at 
fj, — l,t = 0.5. The density wave is absent at ^ = 0, as shown 
in the lower inset. 
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(x(i),3) = 1; ef^ = -1, if mod {x{i),3) = 2; 

e- ' = 0, otherwise. S2/N is the square of the density dif- 
ference between the A and B sublattices. We also define 
the third quasi-structure factor 
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mod (a;(i),3) — 2; e,-"''' = —1, otherwise. S3 measures 
the difference of the boson densities on the B and C sub- 
lattices. In the present model, we find that S2 is always 
zero, which means that the symmetry between the A and 
B sublattices is always kept, thus S3 is also the density 
difference of the bosons on the A and C sublattices. 

The superfluid stiffness is measured in terms of winding 
number fluctuations^'', 



ries expansion (SSE) quantum Monte Carlo method with 
directed loop updatcs^^ for the hardcore bosons on the 
dual of the bowtie lattice. 
quantum phase diagram 

We start with the simple limit = 0. The phase 
diagram is shown in Fig. |31 We find that the system is 
in a phase with striped order II and nonzero superfluid 
stiffness between a Mott insulating phase (MI) and an 
empty phase, except for fi = where the model has an 
exact particle-hole symmetry. 

The phase diagram for 7^ is much richer, see Fig.|4l 
We find SCDW I with p = 1/2 and SCDW II with p = 
2/3 when the hopping is weak, as expected. The phase 
transition between the two SCDW phases is first order. 
A stable SSF emerges outside the two SCDW phases, 
except for p = 0. We will now proceed to discuss the 
phase diagrams in more detail. 

To characterize density wave order, one usually mea- 
sures the static structure factor S'(Q) with Q the 
wavevector 



5(Q)/7V 
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where fc labels sites, and iV is the total number of sites. 
However, in our case, it is more convenient to distinguish 
the two striped phases by introducing the following quasi- 
structure factors, which measure the density differences 
between different sublattices, or, in the language of the 
XXZ model, the various magnetic orders in the z direc- 
tion: 



n^f). 



(4) 



For a = 1, e'l^ = (— 1)"^^*^, where x(i) is the x coordi- 
nate of site i. Si/N measures the square of the density 
difference between the two sublattices consisting of sites 
with even or odd x coordinate, showing a density wave 

(2) 

at wave vector (7r/3,0) if not zero. For a = 2, = 1, 
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where a labels the x or y-direction, and /3 is the inverse 
temperature. Typically the superfluid density ps is the 
average of the two stiffnesses. Considering the anisotropy 
between the x and y-directions, we shall distinguish the 
superfluid stiffness in the two directions. 

We measure the particle density, the quasi-structure 
factors and the superfluid stiffnesses as functions of the 
chemical potential p, and the hopping amplitude t (in 
units of the nearest-neighbor repulsion V , if ^ 0). In 
the simulations, we set ^ Ly ^ L, which is restricted 
to multiples of 6. The total number of sites is = 
L^/2. The inverse temperature was chosen as [3 = 2L/t 
to make sure that the simulations access the ground state 
P 
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FIG. 4. (Color online). Phase diagram with nearest-neighbor 
repulsion present. Blue (dashed) lines indicate continuous 
phase transitions, and the red (bold solid) line represents first 
order transitions. The inset shows the density wave in the 
SCDW I phase at pjV = l,t/V = 0.3, indicated by the filled 
uptriangle. 
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quantum phase diagram for V = 0. 

We first consider the non-repulsive limit = 0. The 
phase diagram is shown in Fig. [31 which is invariant under 
the interchange of particles with holes p — > 1 — p and the 
change of sign /i — ^ —fi. The phase boundaries can be 
well predicted in the single particle picture. The system 
is empty when /i < —lOt/3, and full when /i > lOt/3, 
where 10/3 is the average coordination number (twice 
the ratio of the total numbers of bonds and sites). Our 
simulations confirm this prediction. 

Fig. [5{a) shows the boson density p as a function of the 
chemical potential with hopping t = 1. We see that the 
density varies continuously from to 1 as ^ changes from 
-10/3 to 10/3. In this region, and p^ become nonzero 
with 5*3 > for p > 0, or S'3 < for p < 0, as shown 
in Fig. ^h) and (c). This means the system displays 
both DLRO and ODLRO. To further confirm this, we 
sample the average particle density at the x-th column 

P(^) = {J2i=i "i(2^)/(-^!//2)), i denotes sites in the x-th 
column. A density wave at wavevector (27r/3,0) is seen 
in the upper inset of Fig. [3] for p = l,i = 0.5. The 
density wave for p = — l,t = 0.5 follows by applying the 
transformation p — > 1 — p. For p = 0, no such wave 
appears, as shown in the lower inset. Thus the system 
indeed has a striped order of type II. 

This result can be understood in the following way: 
At p = lOt/S, holes can appear in the system. Each 
hole costs a potential energy p, and gains kinetic energy 
— lOi/3 by hopping freely. The wave functions of the 
holes spread out over the entire lattice, but the proba- 
bility of finding a hole in the C sublattice is 4/3 times 
of that in the B (or A) sublattice due to the coordina- 
tion number difference. The striped density wave thus 
starts to develop. As p decreases, the price of creating a 
hole becomes cheaper and the density of holes becomes 
larger, thus the density wave is getting stronger. As the 
hole density further increases, the hardcore nature of the 
bosons prevents further increase of the density difference 
between the C and A (or B) sublattices. The difference 
approaches a maximum, and then starts to decrease. Fi- 
nally, at p = 0, the lattice is half filled and the difference 
disappears. However, as mentioned above, this striped 
order II does not break the translational symmetry of the 
original lattice, nor the symmetry between the A and B 
sublattices. Thus, we would call this phase bearing both 
ODLRO and DLRO a striped superfiuid (SSF), instead 
of a supcrsolid phase (SS). This SSF phase is a result 
of the interplay of the chemical potential and the lattice 
inhomogeneity and anisotropy. 

This picture can also be understood in the language 
of the XXZ model. At p = lOi/3, the xy components of 
spins start to align, forming a ferromagnetic phase due 
to the competition between the in-plane coupling and the 
external magnetic field h — fi. As the field p decreases, 
the ferromagnetic order in the xy plane gets stronger, and 
meanwhile the magnetization in the z direction decreases. 
The inhomogeneity of the lattice structure makes the lo- 
cal magnetization of z component in C sublattice weaker 
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FIG. 5. (Color online). The particle density p, quasi-structure 
factor 53, and superfiuid stiffnesses (red (light) symbols) 
and pi (blue (deep) symbols) as functions of p at V = 0, t = 1. 



than that in the A and B sublattices. This results in a 
staggered long-range ferromagnetic order in the z direc- 
tion. When the external field p reaches 0, the magneti- 
zation in the z direction disappears completely. 

An interesting phenomenon in this model is that the 
superfiuid stiffness along the a;-direction (perpendicular 
to the stripes) becomes larger than that along the y- 
direction (parallel to the stripes) near half filling. In 
contrast, a larger superfiuid stiffness along the stripes 
than transverse to the stripes was reported in various 
striped supersolid states on the square and triangu- 
lar lattices^iiii^"— , where the lattice itself is isotropic. 
Clearly, the anisotropy of the superfiuid stiffness refiects 
the inequivalence of the x and y-directions of the dual of 
the bowtie lattice, but it is not clear a priori that which 
of them should be larger. 

quantum phase diagram for V ^ 0. 

With the nearest-neighbor repulsion V present, we 
have a different phase diagram, without the symme- 
try when interchanging particles with holes. In the 
single-boson picture, the system becomes empty when 
p < — 10</3. Our simulation results for F > are in 
agreement with this picture. Considering a single freely 
hopping hole on the lattice, one can show that the system 
sits in the p = 1 state when p > \0{t-\-V) which is the 
Mott insulating state. This is true for t/V 3> 1, where 
the single hole added is almost free with its wavefunction 
spreading over the whole lattice. Our simulations show 
that the MI boundary is straight with the expected slope 
in this region. However, the boundary is curved when 
the system is approaching the classical transition point 
p/y = 4,t = 0. For small t/V, holes are created on the C 
sublattice, which costs the same chemical potential p but 
gains a potential energy 4^, more than for those created 
on the A or B sublattices. These holes are forced to hop 
along y-direction. The kinetic energy gain of a single hole 
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FIG. 6. (Color online). The particle density p, quasi-structure 
factors Si/N,S3, and superfluid stiffnesses (red (light) 
symbols) and p| (blue (deep) symbols) as functions of fJ^/V , 
at the cut t/V = 0.3. 



by the second order hopping processes is —At'^/V. Oth- 
erwise, the kinetic energy gain t cannot compensate the 
cost of the potential energy. This explains the stronger 
superfluid stiffness along the y-direction, see Fig. [BJd), 
and results in the curved boundary: /i — AV + /V . 

We now show results obtained by scanning the chem- 
ical potential at constant t/V. In Figs. [6lja)-(d), we 
plot the particle density p, the quasi-structure factors 

/N, 5*3 and the superfluid stiffnesses , as functions 
of the chemical potential ^/V for t/V = 0.3 respectively. 
S2/N is always zero, which means that the boson densi- 
ties in A and B sublattices are always equal. 

Between n/V = 2 and 3.7, the p = 2/3 SCDW II phase 
is found with superfluid density = 0, which means 
that the bosons are localized. The quasi-structure factor 
^3 K, 0.28, which is less than the value 1/3 for the exact 
static SCDW II phase due to the presence of hopping t. 

At p/V = 3.7, the density starts to grow continuously 
as /i increases, indicating a second order phase transition. 
53, and p| are all finite in the region 3.7 < p/V < 4.6. 
The model displays both DLRO and ODLRO. Again, 
since the striped order II does not break the transla- 
tional symmetry of the original lattice, the model turns 
out to be in a SSF phase, which can be understood in 
terms of doping bosons in the SCDW II state. The doped 
bosons are mobile along the y-direction in the C sublat- 
tice and generate stronger superfluid stiffness along the 
y-direction than along the x-direction. The doped bosons 
can not form a domain wall which breaks the striped or- 
der II. Thus the SSF phase is stable. Equivalently, the 
phase can be thought as doping holes in the Mott in- 
sulating state, which hop along the y-direction to avoid 
potential energy cost and form a stronger superfluid stiff- 
ness in y-direction, as described in previous text. 
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FIG. 7. (Color online). The p = 1/2 solid doped with holes, 
(a) Holes (green circles with dashed boundary) added in the 
solid, (b) Lining the holes costs no additional energy, (c) 
A domain wall (dashed red line) is introduced at no cost by 
shifting the right half of the lattice, (d) The holes can hop 
freely across the domain wall and gain kinetic energy. 



A SCDW I phase is clearly observed in the region 
0.68 < fi/V < 2, with p 1/2, Si/N = 0.193, 
5*3 — 0.014, and p'^,p^ converging to 0. Here, the striped 
order I shown in Fig. [T] (a) is slightly adapted by an ad- 
ditional order II. We show the density per row p(x) in 
the inset of Fig. S) The ordering at wavevector (7r/3,0), 
which breaks the translational symmetry of the original 
lattice, is clearly seen. Thus the bosons here form a true 
symmetry-broken solid state. 

Doping this solid (SCDW I) with holes leads to a phase 
separation, instead of supersolidity. A density jump 
caused by a first order phase transition is observed at 
p/V = 0.68, leaving SCDW I for a p < 1/2 superfluid 
phase. The mechanism is similar to what was found for 
the square lattice and the triangular lattice^i^. Adding 
L/2 holes onto the SCDW I solid decreases the density 
infinitesimally in the thermodynamic limit. Each hole 
costs a chemical potential p and gains a kinetic energy 
which is quadratic in t by the second-order hopping pro- 
cesses (see Fig. El^a)). Placing these holes along a verti- 
cal line costs no additional energy, as shown in Fig.lTJb). 
We can shift one half of the lattice to the left (or right) 
by one unit of x coordinate, introducing a domain wall 
which breaks the SCDW I order without cost of energy, 
as shown in Fig. [71[c). By hopping freely across the do- 
main wall, each additional hole gains a kinetic energy —t, 
which lowers the energy of the domain wall state com- 
pared to the bulk supersolid (see Fig. EJd)). Thus, the 
supersolid phase with SCDW I solid order is unstable. 
However, the striped order II due to the inhomogeneity 
and the anisotropy of the dual of the bowtie lattice is 
still present when the boson density p is lower than 1/2. 
We see that 6*3, which shows the density difference be- 
tween B (or A) and C sublattices, is not zero (Fig. IH^c)). 
Since the superfluid density is also nonzero in this region, 
the system is actually in a SSF phase, which finally ends 
when the lattice becomes empty at p/V = 10{t/V)/3, 
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FIG. 8. (Color online). Energy level crossing (a) and multi- 
loop hysteresis (b) for hardcore bosons on a 12 x 12 lattice at 
tjV — 0.2. Three phases with different densities coexist near 
fj,/V = 2. 




FIG. 9. (Color online). The particle density p, quasi-structure 
factors Si/N,S3, and superfluid stiffnesses (red (light) 
symbols) and p| (blue (deep) symbols) as functions of fJ^/V, 
at the cut t/V = 0.4. 



with the line /i = as an exception. Thus the phase 
separation at p = 1/2 is between a SCDW I phase and a 
SSF phase. 

Moreover, we notice that pj is much stronger than p| 
when p is close to 0.5, where S3 ~ 0.05, which is much 
larger than its maximum value at V = 0. This means 
that bosons are partly localized due to the presence of 
a repulsion V. The localized boson density on the A 
(or B) sublattice is larger than that on the C sublattice. 
Holes generated in the A (or B) sublattice contribute 
to the superfluid stiffness along the x-direction, which is 
transverse to the stripes. This anisotropy is again very 
different from what found in various striped supersolid 
states near half filling'* i^i^"—. 

The phase transition between SCDW I and SCDW II is 
rather special. As discussed above, at the transition point 
fi/V = 2, there are actually L/6 + 1 coexisting phases in 
the classical limit. With hopping t nonzero, but small, 
we can still see that these states are degenerate near 
fi/V — 2, as shown in Fig. [8] (a). Away from, but close to, 
fi/V ~ 2, one of the states becomes the ground state, the 
others are metastable. This situation is illustrated by the 
interesting multi-loop hysteresis curves for a 12 x 12 sys- 
tem at t/V = 0.2, see Fig.[5](b). During the simulations, 
we store the last configuration of a finished simulation 
as the initial configuration of the next simulation with a 
new value of the chemical potentiaU^i^. Starting from 
fi/V much less than 2 to ensure that the system stays 
in the ground state, i.e. p = 1/2 SCDW I, we increase 
fi/V and sample the particle density p. The system does 
not jump to the real ground state (p = 2/3) immediately 
when fi/V passes the transition point fi/V — 2. This 
curve is labeled as 'upl' in Fig. IH^b). Then, we start 
with fi/V ^ 2 to ensure that the system stays in the 2/3 



filling SCDW II state, then decrease fi/V. The system 
may jump to the p = 1/2 state which is the ground state, 
closing a hysteresis loop (the curve 'downl'), or jump to a 
metastable phase with 7/12 filling by chance (the curve 
'down2'), when fi/V is small enough. We can use the 
latter configuration as the initial configuration, and in- 
crease fx/V again. It is seen that the system stays in the 
metastable state until fi/V reaches a value much larger 
than 2, as shown by the curve 'up2'. Curves 'down2' and 
'up2' form another hysteresis loop. 

This interesting phenomenon makes the simulations 
near fi/V = 2 very difficult. The data near this point, 
shown in Fig. [51 are obtained by initializing configura- 
tions in the 'right' way. The transition points between 
the SCDW I and SCDW II phases can be found from the 
energy level crossings. 

Increasing hopping t/V larger than 0.37, we find that 
the SSF phase emerges in the region between the two 
SCDW phases. To demonstrate this, we take t/V = 0.4, 
and scan the chemical potential fi. The particle density 
p, quasi-structure factors Si/NjS^, and superfluid stiff- 
ness pf and p| are plotted in Fig. [nia)-(d) respectively. 
Between the SCDW I and SCDW II phases, we see that 
striped order II and superfluidity coexist. It is also clear 
that the transition from SCDW I to the SSF phase is of 
first order, and the transition from SCDW II to the SSF 
phase is continuous. In this region, the superfiuid stiff- 
nesses in the two directions become different with pj > p| 
near half filling. This is also the case for the SSF state 
with p < 1/2 near half filling. For the SSF state with 
2/3<p<l, we see the same anisotropic behavior as 
that at the cut t/V = 0.3. 

The general picture of the phase diagram along fi/V 
for a cut at constant t/V changes a.s t/V becomes even 
larger. At t/V = 0.5 shown in Fig. [TDl the 2/3 fiUing 
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FIG. 10. (Color online). The particle density p, quasi- 
structure factors S\/N,S-i, and superfluid stiffnesses (red 
(light) symbols) and p| (blue (deep) symbols) as functions of 
p,/V, at the cut t/V = 0.5. 



SCDW II phase is the only density wave phase. Although 
the finite size data oi Si/N are nonzero in a large region, 
a finite-size scaling analysis indicates that they finally 
converge to zero as system size turns to infinity. The 
strong hopping destroys the SCDW I order. The SSF 
state with both DLRO and ODLRO is found outside the 
SCDW II phase for ^JL/V < 2.65 and fi/V > 3.5. Close to, 
but outside the SCDW II phase, it is seen that the type 
II striped order persists. Meanwhile, the superfluid stiff- 
nesses start to increase. Hence, the SSF state appears. 
We still see a large anisotropy of superfluid stiffness in 
the SSF states: > p| near half filling, and > for 
2/3< p< 1. 



conclusion 

We have investigated the ground state behavior of 
hardcore bosons on the dual of the bowtie lattice. A 
special state (SSF) with a striped order II (see Fig. [T]) 
and an ODLRO is found even in the absence of repulsion. 
The emergence of such a state is due to the inhomogene- 
ity and the anisotropy of the lattice structure. However, 
neither the translational symmetry nor the symmetry be- 
tween the A and B sublattices of the original lattice is 
broken. Such SSF states should exist on anisotropic and 
inhomogeneous lattice with sites having different coordi- 
nation numbers. Including a nearest-neighbor repulsion 
causes a much richer phase diagram. Two SCDW phases 
with different striped order are found. A SSF phase ex- 
ists outside the two SCDW phases, between the Mott 
insulating phase and the empty phase. Comparing with 
the SSF state with = 0, we see that the striped order 
II is enhanced greatly by the repulsion in the SSF state. 
The phase transition between the two SCDW phases is 
first order. The transition between the SCDW I and 
SSF state is also first order, while the transition between 
SCDW II and SSF state is continuous. We also report 
anisotropics of superfiuid stiffness in the SSF states with 
and without repulsion. In the SSF state with repulsion, 
the superfluid stiffness shows different anisotropics in the 
region near half filling and above 2/3-filling. Hardcore 
bosons on inhomogeneous lattice present rich and inter- 
esting properties. Further theoretical or experimental 
studies are worthy to explore the nature of the SSF state 
in more detail. 
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